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Abstract 

Four-dimensional quantum electrodynamics has been formulated on a hyper- 
cubic Minkowski finite-element lattice. The equations of motion have been 
derived so as to preserve lattice gauge invariance and have been shown to be 
unitary. In addition, species doubling is avoided due to the nonlocality of 
the interactions. The model is used to investigate the lattice current algebra. 
Regularization of the current is shown to arise in a natural and nonarbitrary 
way. The commutators of the lattice current are calculated and shown to 
have the expected qualitative behavior. These lattice results are compared to 

various continuum calculations. 
11.15.Ha,11.15.Tk,12.20.Ds,11.40-q 



Typeset using REVT^ 



* E-mail : dfmiller ©phyast . nhn . uoknor . edu 



1 



I. INTRODUCTION 



The linear finite-element approach has previously been successfully applied to quantum 
electrodynamics The equations of motion comprise a self-consistent lattice gauge the- 
ory, that is, they are gauge covariant to all orders of the lattice spacing. Also, the lattice 
Dirac equation has been shown to be unitary, i.e. the fermion canonical anticommutation 
relations are preserved in time. This is true even in the presence of background electromag- 
netic fields. The lattice Maxwell's equations are shown to preserve unitarity in the absence 
of interactions. The interacting case is less clear due to anomalous behavior which this arti- 
cle will begin to address. In addition, it has been shown that this model avoids the common 
lattice problem of fermion species doubling [||. The no-go theorem 0] is avoided by the 
nonlocality of the interaction term in the Dirac equation. 

The finite-element approach has been successfully applied to massless two-dimensional 
quantum electrodynamics (the Schwinger model) . The axial- vector divergence anomaly 
was calculated to be 



Here, is the axial- vector current, e is the electric charge, M is the number of spatial 
lattice sites, and E is the electric field strength. The derivative is taken according to the 
finite-element prescription: a forward difference in the direction of the derivative, and a 
forward average in the other directions. The relative error is of order as is expected 
from a linear finite-element approach. The anomaly in this model has also recently been 
evaluated by a lattice loop calculation of the vacuum polarization 

In this paper, we calculate the commutators of the vector current in four dimensions. 
These are shown to be roughly consistent with results obtained in the continuum. In addi- 
tion, the regularization of the fermion bilinear ip^ip is shown to be a direct result of requiring 
gauge invariance and unitarity. 

In Sec. |I| we briefly review several continuum calculations of these commutators. The 



ambiguity inherent in these results is illustrated. In Sec. Ill our model is introduced in the 



form of the four-dimensional lattice equations of motion. Several important features of these 
are presented. In Sec. we present results for vector-vector commutators, and in Sec. [V| 
we present our conclusions. 



II. CURRENT COMMUTATOR ANOMALIES IN THE CONTINUUM 



It is well known that there are anomalies in the commutators of the electromagnetic 
currents. Schwinger first calculated a non-zero value for the commutator of the charge 
density and vector current using a point-splitting regularization of the current [0, 



[X) = limeip^ (x — -e) 7°7^'?/^ (x + -e 



(2.1) 



where e is space-like and the limit is taken symmetrically in space. Using this current, he 
arrives at a value for the commutator of 



J° (0,x), J(O)l) = -SV5 (x) 



(2.2) 



where S is the divergent limit of (2e^/37r^)(e^)^^. 

Later calculations [^|-[Tl[] revealed an additional finite contribution to this commutator. 
Its value depends on the type of regularization chosen to regulate the bilinear current. Hence, 
Eq. ( p.2|) is modified to read 



^( f J° (0, x) , J (0)1 ) = -SV6 (x) - 4 VV2(5 (x) 

L J 71"^ 



(2.3) 



where S is the divergent Schwinger term and d depends on the method of calculation. 

The Bjorken- Johnson-Low method defines the equal-time commutator in terms of the 
high-energy limit of the time-ordered product 



lim -ipo I d^xe''"'{TAix)Bm 

O^oo J 

= J d^xe-'P^'HA (0,x),5(0)]). 



(2.4) 



This is shown PJTT|] to be equivalent to the following position-space prescription: 
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{[A (0,x),B(0)]) 

= lim [{A (r^, x) B (0)) - {B (0) A {-r^, x))] . (2.5) 

Using this prescription, the value of d in Eq. { \i.'d\) is calculated to be = 1/12 [|T^,|ri[|. 

The value of the finite piece of Eq. ( |2.3| ) can also be calculated using a point-split current 
such as defined in Eq. ( p.l|) . Boulware and Jackiw use this to derive a value of d = 1/60 
||To| . These authors also use a generalized point-split current, allowing e to be a general four- 
vector, to derive a value of d = 1/96. To further illustrate the ambiguity in this calculation, 
we give the results of Chanowitz [|ll|]. He shows that is one uses an unsymmetrical definition 
of the current 

J^" (x) = lim {x + e) 7 VV^ (x) , (2.6) 

with e spacelike, then a value of = 1/15 is derived. 

Evidently, the results for Eq. ( p.3|) depend on the definition one chooses for the regulation 
of the current. Of course, this is to be expected, since the quantity is divergent. Presumably, 
the exact value of the coefficients of the Schwinger terms will have no physical consequence. 
However, their very existence proves to be problematic in establishing the formal properties 
of the theory. 



III. THE LATTICE EQUATIONS OF MOTION AND REGULARIZATION OF 
THE ELECTROMAGNETIC CURRENT 

A. Lattice Dirac equation 



The finite-element formulation of the gauge-invariant Dirac equation was carried out 
in Refs. ||l|J^. We will merely quote the results here, expressed in terms of the spatially 
averaged fermionic field: 



E- E 

I m'.=l m'=m,+l 



(-1) 



m±,n 
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n 



2^7- 



j M 



X! -^mjL\rnjrn'.,n'^m'.,m^,n — 0, (3.1) 



n'=0 



m . =1 



where a sum over the repeated index j is understood, and the overbar represents a forward 
average over that coordinate: 



(3^2) 



(Recall [m that with M odd, is periodic on the spatial lattice.) We have chosen a hypercubic 
lattice with lattice spacing A. 

Eq. ( p.l|) was explicitly constructed to be invariant under the local gauge transformation 



40 _|_ _ { \— _ A_ 



(3.3) 



This particular choice of transformation is made so that the mass term in the Dirac equation 
( ^yil) is automatically covariant. 

To simplify things, we choose the temporal gauge (A° = 0). In this gauge, X° vanishes 
identically and Eq. ( p.l|) contains only fields at time n. The explicit form of the spatial 
parts of the interaction is 



mx ;nij ,m^ ,n rnj,mj \ j 



M 



1 + COS sgn (m" - mA sgn \m\ - m'A Cx,m",n sec C^^^„ 



, m"=l 



+ « (—1) ^ ^ sm 



/ M 

\ 



IVl I 

J2 sgn (m^' - m^) sgn (m'^ - m'j) Cx,m",n sec C4x,n- (3-4) 



my=i 



We have used the abbreviations 



eA 



/-J — A 

Smj_,mj,n <-> m,,m„-l,n' 



Tn±,'n 



M 



(3.5) 



and 
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sgn [X) 



+1, X > 
-1, x<0 

+ 1, rrij > m'- 



0, rrij 



m 



■1, rrij < m'j. 



Eq. ( |3.1| ) can be written in a more compact matrix notation: 



A 



(Q - In) (j)n + (^(Pn = 0. 



Here, Q and I are matrices in the spatial indices 



.In) , ^mx,m' I ;m,,m'.,n' 



and 0n is a vector of the spatially-averaged fields 

Notice that Q is related to the lattice derivative, 
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and the lattice version of the covariant derivative is 



I (Q - I„) ^ T>. 



Eq. (|3.8| ) can be solved for the transfer matrix defined by 
giving 



_ l-7°7-(Q-In)+^^7° 
" l + 7°7-(Q-In)-^^7°' ^ ^ 

with V = From Eqs. (|3.4|) and (|3.9|) , it is clear that both Q and I„ are anti-hermitian. 
Therefore, the transfer matrix defined in Eq. ( p.l4| ) is unitary. Hence, the spatially averaged 
fields (pn can be taken to be the canonical fermionic fields. 

B. Lattice Maxwell equations 



Using the notation of Eq. ( p.9|) , the lattice Maxwell's equations are written as 

— Q • = J°, 

E„_|_l = J„ — Q ■ n, 

:F^ = -^{Q'AL-Q^Al), (3.15) 



where 

E = E— 

An = AnYn, (3.16) 

Gauge invariance of Eqs. ( p.l5| ) under is assured as long as the current J^^^ is con- 
structed to be gauge invariant. 

C. Current regularization 

The simplest choice of gauge-invariant current is 

J^,n = e?M,7i7''fe (3.17) 

= e0ni.n7''0m,n, (3.18) 



which is invariant under the transformation ( p.3[ ). Notice that the covariant field ^jjm,n that 
is involved in the definition of is not the same as the canonical field 0m,n defined in 
Eq. ( |3.13|) — they differ by a temporal averaging. This forces point-splitting of the current: 



m',m" 

(3.19) 

The important point is that this is not introduced in an arbitrary way. The regularization 
of the current is mandated by requiring gauge-invariance and unitarity of the equations of 
motion. 



\ / m'.m \ / m.m'' 



IV. VECTOR- VECTOR COMMUTATORS ON THE LATTICE 



A. Analytical results 



Using the definition of the electromagnetic current Eq. ( p.l8 ), the commutators of various 



components can be calculated. To lowest order in e, the commutators can be calculated using 
the zeroth-order contribution of the transfer matrix 

_ l-7°(7-Q) + ^^7° ...^ 
" 1 + 70 (7 ■ Q) - iv^^ ^ ■ ' 

and a free field Fock space expansion for the Dirac fields 

"yM5('^)g27rjp-m/M 



(T,p 




P P 



where the spinors are normalized according to 



^,M,(.V^^^^7!^i(7it)+^ (4.3) 



with 



tp 



tan 



M 



UJ-r 




(4.4) 
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Then the canonical anticommutation relations for the Dirac fields 



are satisfied if 



{</'m,n,0L',n} = ^^m,m' (4.5) 



= — ^(5pp,5"'"' 
{^P'<'} = (]^W^^'^', (4-6) 



and all other anticommutators of these operators vanish. 
The results of the calculations are presented below: 



tO tO 



(4.7) 



JLn, JL'n ) = (4.8) 



47>2 tan [ ^ 

7O jj \ ^ V V / 27rip (m-m')/M ( a n\ 



B. Numerical evaluation 

The matrix (in m and m') represented by Eq. ( [4.9| ) is shown for lattice size M = 
9 and mass fi = in Fig. |l[ The abscissas of the plot correspond to one-dimensional 
representations of the three components of m and m', respectively. (That is, the base-M 
number (mi, m2, ^3) is converted to a base-10 value of the abscissa.) It is apparent that the 
leading-order behavior of this commutator is a first derivative of the Dirac delta function, 
as expected. 

To make this comparison more quantitative, we fit this result to the functional form 
expected in the continuum, Eq. ( p.3| ). The lattice analog of the Dirac delta function is 



6 (x) ^ -J-(5^o = -^-^ y e^-'P-™/*^ (4.10) 
^ ' A3 (MA)' p 



so we take as a trial function 



tO jj 



2r 27rip-(m-m')/Af 



^ y p. y a2.+lP. (A) (-1)^ ^ 2 

(MA)'MA^ r=o ^^^^ 

(4.11) 

The coefficients a do not depend on the lattice spacing A (for mass /i = 0); they presumably 
remain finite in the continuum limit, A — 0. The continuum behavior of the terms in ( [4.11 ) 



is dictated by the functions Pr (A). These are inserted to force each term in the series to be 
of the same order in A as the commutator ([4.91), namely 1/A®: 



P,(A)=A'"-^ (4.12) 

Therefore, in the continuum limit, the ffist term in Eq. ( [4.11[ ) will be quadratically divergent, 
the second will be finite, and the rest will vanish, as expected. 

To perform the fit, the data is ffist converted from a three-dimensional to a two- 
dimensional representation. (This is necessary because the generation of 3-D plots is im- 
practical beyond M ~ 9 and because 3-D ffis are prone to difficulties.) Specifically, a plot 
such as Fig. |1| is projected onto a plane orthogonal to the m, m' plane and to the diagonal 
of the matrix. This yields a graph such as the solid curve in Fig. 0(a). The curve is then ffi 
to a similar projection of the trial function in Eq. ( [4. lip . In particular. Fig. ^ shows fit vs. 



data for a ffi of the commutator (^]9[) to a linear superposition of three spectra, i.e., R = 3 



m 



It turns out that a ffi to five spectra produces comparable results. The addition of the 
last two spectra neither particularly improves the ffi nor radically changes the coefficients 
of the original three. (That is, the ffi is stable.) In addition, the coefficients of all but 
the ffist two spectra are vanishingly small. In Fig. ^ we show the results for a ffi to ([4.11 



with i? = 3, 4, and 5 for various lattice sizes. These are compared to some of the results of 
continuum calculations. The lattice values, though of the same order as other results, are 
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significantly different. This is as expected, since (|3.19|) is yet another definition of current 



regularization. (Note that the result of Chanowitz |TT| is within one sigma of our result. 
This is interesting in that both his method and ours use a temporal averaging to regularize 
the current.) 

A note should be added to discuss the parameterization of errors in the data: Since 
the data (|]9|) is computer-generated, any "measurement" errors should be strictly due to 
round-off errors. Presumably this will be roughly the same for each data point. Therefore, 
it is reasonable to assume that all measurement errors are equal, = cr^. If we further 
assume normally distributed errors, then the value of for the fit is 

X^ = ^ (data-fit)^ (4.13) 

data points 

If we now require the fit to be "good" , i.e x'^ per degree of freedom equal to one, then we 
can determine a value for cx^, 

^2 _ Sdata points (data - fit) ^ 

degrees of freedom ' 

where the number of degrees of freedom is the number of data points less the number of terms 
in the trial function. Multiplying this by the diagonal elements of the covariance matrix C 
(which is determined as a by-product of our minimization process), gives the variances of 
the coefficients 

(ai) = a^Cjj. (4.15) 

These are plotted along with the values of the coefficients in Fig. 

There are a few features of the graphs of Fig. |^ which merit some discussion. First, there 
is the consistent kink in the data at lattice size M = 13. We have no explanation for this 
behavior except to suppose it is just a numerical quirk of our code. 

Second, the scattering of the data points increases significantly with the number of fit 
spectra R. This seems to be a problem with numerical precision in the fitting procedure — 
Each successive term in the trial function ([4.11|) is an order of magnitude or two larger than 
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the previous one. Thus, at higher numbers of fit spectra (i? ~ 5), the range in values of 
terms in ( |4.11| ) approaches the numerical precision of the data. For example, if we continue 
to R = Q, our procedure breaks down — the coefficients generated by the fitting algorithm 
do not give a good fit to the data. 

Lastly, the errors associated with the data do not seem to account for the scatter of 
the points (especially for _R = 5.) Remember, though, that the derivation used to give the 
errors in our data is based upon several assumptions. These are not rigorously justified. 
Rather these suppositions serve to compensate for our lack of knowledge of the fundamental 
source of errors. In order to summarize the data, incorporating the errors in a more model- 
independent manner, we have averaged the data points at all lattice sizes (excluding M = 13) 
and plotted these results alongside the other data. (Figure |^) 

The data for i? = 3 and R = 4 are in good agreement over a wide range of lattice sizes. 
These only disagree at low M, where the model is not expected to yield good results, and at 
M ~ 31, where the fitting procedure for R = 4 begins to break down. This agreement gives 
support to our claim that the fit is stable with respect to the number of fit spectra R. Thus, 
we are able to consider the best data (with i? = 3) as representative of the commutator 
( [4.91 ). The results in the fit of the trial function ( [4.11| ) with i? = 3 to the commutator (^ 
are given in Fig. ^. 



V. CONCLUSIONS 

Using the definition of the current (|3.19|) , the current commutators are calculated and 
are seen to be anomalous. These exhibit the qualitative behavior expected from continuum 
calculations, having both quadratically divergent and finite contributions. The coefficients 
of these are of the same order of magnitude as the various continuum calculations. In fact, 
our results seem to be in quantitative agreement with those of Chanowitz [|ri|] who also uses 
a temporal point-splitting regularization scheme. 

Furthermore, the results arise from a regularization of the current that is both natural 
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and nonarbitrary. The requirements of gauge invariance and unitarity dictate a current 
which exhibits a temporal averaging of the composite fermionic fields. This averaging is 
distinct, yet akin, to a point-splitting regularization. 

We have also begun work on vector-axial- vector commutators. These are more interesting 
since they can be directly related to the chiral anomaly |]T^ 

^(^) = 4^e°^^'{^„F,,}, (5.1) 

by the Gauss' Law constraint, 

V ■ E = J°. (5.2) 
Preliminary results are promising, but much work remains to be done. 
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FIGURES 

FIG. 1. Three-dimensional plot of for M = 9 and // = 0. 
FIG. 2. (a) Fit vs. data of two-dimensional projection of (4.9) and (4.11) with R = 5. {M = 21 



and ^ = 0) (b) Smaller scale picture of the same. (The position variable has been normalized to 
range from —1 to -|-1.) 

FIG. 3. Coefficients of a 3, 4, and 5 spectra fit for different lattice sizes. (// = 0) (a) Coefficient 
of first derivative of delta function, (b) Coefficient of third derivative of delta functions. The 
right-most data points summarize the data excluding M = 13. Also plotted are various continuum 
results for d in Eq. (|2.3|). 



FIG. 4. Summary of the results for a three spectra fit. 
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